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Abstract 

Self-consistent, mean-field description of charge injection into a dielectric medium is modified to 
account for discreteness of charge carriers. The improved scheme includes both the Schottky barrier 
lowering due to the individual image charge and the barrier change due to the field penetration 
into the injecting electrode that ensures validity of the model at both high and low injection 
rates including the barrier dominated and the space-charge dominated regimes. Comparison of 
the theory with experiment on a unipolar indium tin oxide/poly(phenylene vinylene)/Au device is 
presented. 

PACS numbers: 72.80.Sk, 72.80.Le, 73.30.+y 
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I. INTRODUCTION 



Injection of charge carriers from a conductor to an insulating medium has a crucial impact 
on the functioning of a wide range of electronic devices: light-emitting diodes, field-effect 
transistors, capacitors etc. using inorganic or organic semiconductors, wide-gap insulators 
and ferroelectrics. Nevertheless, the description of this process still remains controversial. 
By the simulation of inorganic and organic semiconductor devices the fixed values of charge 
carrier densities at the boundary with the insulating region are often used as boundary 
conditions^-. For charge transport in bulk dielectrics the electric field at the injecting inter- 
face is often taken equal to zero assuming space-charge limitation of the current^ - — . Given 
finite values of the electric field at the interface are supposed in numerous papers simulat- 
ing injection as Fowler-Nordheim (FN) tunneling through the surface energetic barrier— ~— 
or Richardson-Schottky (RS) thermionic emission over this barrier—^"—. However, the 
problem is that the field and the carrier density at the interface are, in fact, unknown. 

Moreover, the injection treatment in the spirit of the FN or the RS approach applies, 
strictly speaking, only to the case of a single particle penetration since the presence of other 
particles immediately changes the electrostatic potential profile used in these models. This 
makes the application of the mentioned models questionable, particularly at high injection 
levels. A consequence of the single-particle approach is that the characteristics of a dielectric 
medium are virtually absent in the results. Expressions for the FN and RS injection currents 
look as if it is injected into vacuum (except for the dielectric permittivity) and do not 
depend either on density of states available or on the charge carrier mobility in the dielectric 
medium, that contradicts experiment^. This deficiency of the RS model was properly 
noticed in Refs.- 1 ^ where the drift-diffusion equation for charge carriers in the insulating 
medium was combined with RS injection to introduce the carrier mobility and the density of 
states in the injection current, however, the question remains whether a boundary condition 
derived from a single-particle model may be coupled with many-particle equations. In view 
of disordered insulators, sophisticated extensions of the RS approach describing injection as 
random hopping of a charge carrier in the same potential profile with account of possible 
injection into the tail states below the barrier were developecr^^i 1 ^ - — which, however, do 
not fix the general deficiency of the model of injection which remains a single-particle one. 

For a proper description of the charge carrier injection a self-consistent determination of 
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the field and the carrier density at the conductor /insulator interface is necessary which would 
contain both limits of weak (single-particle) and strong (many-particle) injection including 
the space-charge limit. Recently, a self-consistent continuous description of injection in terms 
of carrier densities and mean fields was advancecr^^ based on the matching of the electric 
displacement and the electrochemical potential at the interface. This one- dimensional treat- 
ment implies averaging of all variables over the plane perpendicular to the injection current 
density which pertains to relatively high carrier densities. The model exhibits a plausible 
crossover from the barrier dominated behavior at low voltages to the space-charge dominated 
behavior at high voltages and also a field-induced reduction of the injection barrier. Never- 
theless, at medium voltages, this reduction is much less than the reduction of the Schottky 
barrier due to the single-particle image potential of three-dimensional nature^ which seems 
to have been confirmed in experiments^ 1 ^ 1 ^. Apparently, there are restrictions on applica- 
tion of the mean-field description of the charge injection at low carrier densities. 

The importance of proper account of the discrete, single-particle nature of the interaction 
between an injected carrier and the injecting electrode was emphasized in the Refs.— ~— 
where a sophisticated numerical treatment of the one- dimensional hopping transport in 
organic semiconductors was developed including the contribution of the modified single- 
particle image charge potential. However, for the boundary conditions in this model, the 
tunneling from an electrode was employed assuming numerous atomistic parameters of the 
involved media instead of the density of states available and the carrier mobility. 

The most advanced approach to charge injection including the Schottky barrier lowering 
and the space charge effect was recently developed in Ref.— where three-dimensional (3D) 
hopping of charge carriers over the sites of a cubic lattice with randomly distributed energy 
levels was considered as well as a sophisticated one-dimensional (ID) continuous model. 
The site occupancies and the electric field were calculated self-consistently by solving the 
three-dimensional master equation and the Poisson equation in consequent iterations with 
account of the field-dependent injection barriers. On the other hand, the contribution of 
the individual image potential into injection was accounted inconsistently which lead to an 
overestimation of the image effect as will be shown below. 

The primary aims of this paper are to establish the requirements, which allow for the 
mean-field treatment of charge injection, and to extend the boundary conditions for the 
continuous description^^^ in such a way that the effect of discreteness of charge carriers 
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is considered. We first look for a simple criterion separating the regions of system parameters 
(voltage and injection barrier) where either the single-particle (SP) or the many-particle, 
mean-field (MF) concepts are valid. Then we formulate self-consistent boundary conditions 
which account for the discrete character of the particle interaction with the electrode and 
provide the crossover from low to high injection regime including the case of space-charge 
domination. Finally, the results are compared with an experiment. 

II. COMPARING A SINGLE-PARTICLE AND A MEAN-FIELD CONCEPTS 

Adequacy of the MF description of the charge carrier injection or, alternatively, of the SP 
picture is determined mostly by the relation between two respective characteristic lengths. 
In the first approach this length is given by an average, three-dimensional distance between 

— 1/3 

the injected charge carriers defined as r s = p s where p s is the density of these carriers. 
It is assumed here that the distribution of the particles is homogeneous at least over the 
distance of few r s in all space directions. In the second approach the relevant length is the 
distance x m from the plane surface of an electrode to the position of the potential maximum 
of the Schottky energetic barrier—. 

We verify now the validity of the SP approach from the point of view of the MF approx- 
imation and vice versa. If the density of injected carriers p s predicted in the continuous 
consideration is so high that r s <C x m then the presence of many particles between the 
surface of the electrode and the maximum of the single-particle potential makes the SP 
approach inappropriate; nota bene that MF model density p s is used in this estimation be- 
cause in the consistent SP calculations the density of injected carriers makes no sense. If, 
on the other hand, x m <C r s then the interaction of a single injected charge carrier with 
its image is much stronger than with the other distant injected carriers which makes the 
MF approximation inappropriate. At first sight, the criterion x m ~ r s must roughly sepa- 
rate the regions of relevance of the discrete SP and the continuous MF concepts^. In fact, 
comparison of the SP image force exerted upon an injected carrier with random forces due 
to other charge carriers needs more sophisticated analysis and involves additionally other 
characteristic lengths. To perform this analysis we first delineate below the SP and the MF 
approaches and then study their marginal validity. To make this delineation quantitative 
we consider in the following exemplarily injection of holes from an indium tin oxide (ITO) 
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electrode into an organic semiconductor. 



A. Boundary conditions in the mean-field approach 



Here the treatment in terms of continuous carrier densities and fields is assumed. In the 
linearized Thomas-Fermi approximation^ the electric field in the electrode, occupying the 
space at x < 0, reads 3 ^- 



F c {x) 



J 



F c (0) - ±- 



O, 



X \ J 

ex P 7 — + — > x < 0, 
Itf J o c 



(1) 



where j is the position-independent steady-state current density, a c the conductivity of the 
electrode, and 
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is the Thomas-Fermi screening length. Here e is the positive elementary charge, eo the 
dielectric permittivity of vacuum, e c the static relative permittivity of the electrode, 
and are, respectively, the equilibrium values of the carrier density and of the chemical 
potential in the electrode far away from the interface^. 
The electrochemical potential in the electrode equals^ 
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where Ej, terms the position of the conduction band bottom and <j> c (x) the electrostatic 
potential in the electrode. 

The electric field in the organic semiconductor, F s (x), obeys equation^^ 
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where fi s denotes the hole mobility in the organic semiconductor, e s its static relative per- 
mittivity, k the Boltzmann constant and T the absolute temperature. Nonlinear equation 
OH) is usually solved numerically^ but can also be solved analytically 33 ' 3 ^ which offers a 
certain advantage as can be seen below. The solution reads 
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where Ai and Bi denote the Airy functions^, and i = je 2 lj> F / (2fi s eoe s (kT) 2 ). For the semi- 
infinite geometry considered in this section the constant Ci should be chosen equal to zero 
since the charge carrier density vanishes asymptotically^ whilst the constant C\ must be 
determined from the boundary condition at the interface. 

Assuming Boltzmann statistics for charge carriers in the organic semiconductor we pro- 
ceed, similar to Sze for inorganic semiconductors^, with the electrochemical potential 



k s (x) — E b + Koo + A + kT In 
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+ e<fi s (x), x > 0, (6) 



where the zero- field barrier A is here given by the difference between Ej, + and the highest 
occupied molecular orbital (HOMO) in the organic medium, N is the density of states for 
holes at the HOMO-level, p s (x) the density of holes and 4> s (x) the electrostatic potential in 
the organic medium. In fact, in disordered (organic) semiconductors the injection barrier 
becomes an ill-defined quantity because such systems do not exhibit sharp band edges. 
The barrier is then often defined as a difference between the Fermi level in the injecting 
electrode and the center of the density of states (DOS) distribution at the HOMO-level 
which is typically of the Gauss type^^^'^i'^'i&^I. Using in the following the narrow-band 
approximation we assume a negligible width of the mentioned Gauss DOS a <C A. We 
note, on the other hand, that the following results apply for wide-band non-degenerate 
semiconductors and insulators as well in which case 2(2nmkT / h 2 ) 3 ^ 2 with m the effective 
mass of carriers, and h the Planck constant should be taken for the effective density of 
states^^. 

To obtain the boundary condition at the interface at x = we assume that 

1. electrostatic potential is continuous, C (— 0) = S (+O), which means absence of a 
dipole layer at the interface; 

2. electric displacement is continuous, e c F c (— 0) = e s F s (+0), which implies absence of a 
surface charge at the interface; 

3. electrochemical potential is continuous through the system, which provides matching 
at the interface, k c (— 0) = k s (+0). 

Note that the requirements 1. and 2. account for the image charges of all injected carriers in 
the mean- field approximation. The conditions 1.-3. together deliver the boundary condition 
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at the interface^ 1 ^ 

which connects the values of the field and the carrier density in the organic semiconductor at 
the interface. Due to presence of the current this boundary condition describes, in principle, 
a non-equilibrium state of the contact but in fact the contribution of the current in Eq. (J7J) 
can be neglected here because it is very small in all practical cases concerning organic 
semiconductors^ 1 ^. That is why it will be omitted in the following modifications of the 
formula ([7]). We note also that the condition ([7]) does not define explicitly the values of the 
field and the carrier density at the interface; these quantities can only be found by satisfying 
Eq. ([7D together with Eq. 

Since Boltzmann statistics is assumed to calculate the chemical potential in the organic 
medium the argument of the exponential function in the formula ([7]) must be negative which 
imposes a restriction on the field value F s (0) < Fi im = e c A/ (e s elTF)- The form of Eq. (J7j) 
implies an effective barrier for hole injection 

A e// = A - eF s {0)l TF e s /e c = A - eF c {0)l TF , (8) 

which states a barrier modification due to the electric field at the interface. The barrier 
correction in Eq. (|H]) is linear in field and may be positive or negative depending on the 
sign of the local field at the interface. The physical reason for this barrier modification is a 
contribution of the work done by the mean electric field on electrons in the electrode. 
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B. Schottky barrier and image potential in the single-particle approach 

In the classical description of a Schottky barrier— i22i22 ) the barrier for the emission of 
a single particle through a plane interface is formed by the superposition of the zero-field 
barrier A, the potential of attraction to the image charge of the opposite sign at the electrode 
and the contribution of the mean electric field in the organic semiconductor F s (x): 

e 2 

U(x) = A- — eF s (0)x, x>0, (9) 

16iTe e s x 

where x denotes the distance from the ITO/organic interface, and F s (0) the value of the 
electric field at the interface. The electric field F s (x) is assumed positive and virtually 
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constant within the distance ~ x m from the interface, since for the negative field the potential 
exhibits no maximum and does not allow particle escape from the potential well. The 
maximum of the potential barrier, Eq. (jUJ), is located at a position 



where the latter provides the so-called Schottky-barrier lowering . This barrier modification 
is proportional to the square root of the field and always reduces the barrier in contrast to 
the barrier correction in Eq. (JS} which is linear in the field and may have different signs. 
The physical reason for the Schottky-barrier lowering is also different: in contrast to the 
contribution of the mean electric field at the side of the electrode as described by Eq. (jSJ), 
Eq. (llljl results from the energy profile modification in the dielectric medium due to the 
individual image-charge effect. 

Eq. (Q implies ideal screening of the electrostatic field by the metal surface at x — 0. 
ITO is, however, a highly doped semiconductor which does not provide ideal screening of the 
field. Taking into account the dielectric permittivity of the electrode e c and the characteristic 
length of the field penetration into the electrode Itf the energy of the injected charge carrier 
is substantially modified 4 ^. While at large distances from the interface x 3> a = (e s /e c )lTF 
the energy approaches asymptotically Eq. (Q, at small distances iCa, the Coulomb term 
~ e 2 (e s — e c )/((e s + e c )167re e s x) prevails, which can even change sign from attraction to 
repulsion depending on the magnitudes of respective permittivities. At such small distances 
a consistent quantum-mechanical treatment of injection becomes necessary which, in turn, 
essentially modifies the energy profile in the close vicinity of the interface^ 1 ^. For our 
consideration, the important maximum of the potential (J9J) would fall into the region x m ~ a 
for a voltage of about eL/167reoe s a 2 with L the device thickness. For typical parameters of 
both involved media given in the Tableland L ~100 nm this voltage amounts to 200 V which 
is far too large for organic materials. Therefore in the following the classical approximation 
for the image potential, Eq. (jnj), will be used keeping in mind possible modifications of this 
interaction when other electrode and dielectric materials are involved. 
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and amounts to U(x m ) = A — e5<psch with 
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C. Microscopic comparison of the single-particle image force and stochastic forces 

To establish on what terms the SP approach transforms to the MF one we first generalize 
the formula (Q. In classical statistical mechanics the microscopic density of injected particles 
is defined as 4 ^ 

p(r,t)=]£<5(r-r,(t)), (12) 

i 

where 5(r) is the three-dimensional Dirac delta- function and vectors rj(£) indicate random 
positions of all injected particles and their images. The energy of an injected particle is then 
given by 

U M (r,t) = U (t) - + / dr'p(r>,t) e2sg * {x,) (13) 

167re e s x J 47re e s |r - r'| 

and the x— component of the microscopic force exerted upon this particle reads 

dU e 2 



dx lQ-weQEsX 2 

A (dvW,t) j^&l (14) 
dx J ' 47reoe s |r — r'l ' 



where the first term in both Eqs. ( {TBI and (1T41) presents the contribution of an external field, 
the second term is due to the individual image charge and the third one accounts for the 
other injected particles. The signum function sgn(x') in the last term accounts for opposite 
sign of the image charges. To judge whether the SP approach is relevant, the impact of 
the deterministic individual image effect should be compared with the impact of the other 
particles. Due to the stochastic nature of the interaction of the individual particles leading 
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to the energy and force fluctuations the comparison must be performed in two ways: the 
mean value or, alternatively, the variance of the third term in Eq. f lT3|) or Eq. (|14p should 
be compared with the second term. As long as the mean value and the variance of the 
stochastic term is much less than the individual image charge contribution the SP approach 
prevails. 

1. Mean-force criterion 

By configurational averaging of Eq. (|14p over all possible charge carrier positions the 
mean, macroscopic density of particles becomes a continuous function of the only variable 
x: < p(r,t) >= p s (x). After integration over the space variables y' and z' the mean force 
results in the form 



where the second and the third terms should be compared at x = x m . Keeping in mind the 
hopping transport over the random atomic or molecular sites typical of disordered (organic) 
semiconductors, we assume here that the characteristic length in the space dependence of 
p s (x) is much larger than the nearest-neighbor distance between the sites occupied by the 
charge carriers, tq ~ 0.1—1 nm. 

Eqs. (fT3|) - (|T5|) make sense if the density of injected particles decreases at distances much 
smaller than the device thickness, L, otherwise the electrostatic potential has to be modified 
to account for the other electrode as it is done in the Appendix. Nevertheless, even over- 
estimation of the integral in Eq. ( 1T51) for the very low constant concentration p s (x) ~ r~ 3 
at — L < x < L with r s ^> L results, by comparison of the second and the third terms, in 
assessment x m ~ a/ r^/8nL ^> L. Since such x m cannot be realized, it means that the third 
term in Eq. f lT5|) remains always much less than the second one, so that the SP contribution 
dominates in the limit r s ^> L. This formal statement is in agreement with the obvious fact 
that, for r s > L, the description in terms of the continuous charge density is not valid. 

For higher concentrations with r s <C L we may apply to the evaluation of the integral 
in Eq. (fT5|) the equilibrium solution of Eq. @ which is known to prevail in a wide range of 
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applied fields 3 ^-^ 



Ps{x) =p s {0)/{l + x/XV2) 2 , 



(16) 



with the Debye length A = y^e e s kT /e 2 p s (0). By comparing the second and the third terms 
in Eq. (115fl the criterion results 



*-=W® ,+r Vi- (17) 



where one more characteristic length = e 2 /327reoe s fcT = r 3 /327rA 2 of the order of the 
Coulomb capture radius 50 appears. For typical parameters of organic semiconductors at 
room temperature 3 *^ 2 - l T is about 2 nm. In the range of concentrations L~ 3 <C p s (0) <C 
it follows from Eq. (ITT)) in a good approximation that 



This means that, for x m smaller than that given by Eq. (I18j) . the SP term in the force, 
Eq. ( !T5|) . dominates while in the opposite case the SP contribution may be considered as 
embedded in the MF term. 

For higher concentrations p s (0) S> l^? the equilibrium density ( TT6|) does not apply any- 
more. By increasing external field the injected charge density is known to transform from 
the diffusion induced equilibrium distribution (TIB"]) to virtually constant distribution over 
the device 3 ^ which transforms by further field increasing to the distribution p s (x) ~ l/\/^ 
typical of space charge limiting currents^. Assuming that at r s < Ij- the concentration 
p s {x) ~ l^ 3 remains constant all over the device and x m <C L one obtains, by equating the 
second and the third term in Eq. ( fl5l) . the criterion x m ~ r s \Jr s /\QirL ~ lO -1 /^- According 
to Eq. ( FlOl) this can only be realized for very high fields where the classical SP approximation 
will be violated because of quantum effects discussed in Sec. IHBI Thus, in the region of 
realistic electric fields the SP contribution remains much lower than the MF one so that the 
MF approximation prevails for r s < It- 

For even higher densities, when r s becomes much less than Z^, the MF approximation of 
Sec. Ill Al is expected to fail because of violation of the criterion of applicability of Boltzmann 
statistics, F s (0) < Fu m . On the other hand, the classical SP approximation might also be 
violated in this region because of the mentioned quantum effects. In any case, this all 
happens at rather high voltages as will be seen in the following. 
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2. Energy fluctuation criterion 



Considering the variation of the microscopic energy (fl3|) we will need a second central 
moment of the microscopic carrier density which can be written in the forrn^ 

< Sp(r',t)5p(r",t) >= p s (x)5(r' - r") +p 2 s (x)g 2 (r' - r"), (19) 

with 5p(r, t) = p(r,t) — p s {x) and g 2 {^) the correlation function. In terms of hopping over 
the molecular sites, this form assumes that p s {x) changes on the scale much larger than 
both the distance between the nearest-neighbor sites, tq, and the characteristic length of the 
correlations contained in the function g 2 {r). In Refs.— ^ an assumption of the short-range 
correlations of the charge carrier positions was adopted which only excludes the double 
occupation of sites. This implies that g 2 (r) = — 1 when |r| < r and g 2 {r) = otherwise as 
in the model of hard spheres^. Such space correlations are compatible with the hypothesis 
of the uncorrelated Gauss disorder of the energy levels at different sites in a random hopping 
system considered in Refs.— ^L^k, The validity of the short-range space correlations of the 
particle positions can be verified by comparison with the screened Coulomb correlations 
given by the function^ 

"■'-l, (20) 



g 2 (r) = exp 



8l T ( r 



which is also appropriate for the charged plasma of injected particles of the same sign. 
This formula defines indeed short-range correlations for medium and high charge densities 
when r s < 10 nm since in this case A < 2nm. However, for lower particle concentrations 
with r s ~ L/2 ~ 50 nm the correlation length becomes about A ~ 20 nm. Other possible 
correlations in organic systems and their effect on the results presented here will be discussed 
later in the concluding section. 

The variation of the energy (fl3|) is due to the stochastic term only and can be estimated 

as 

< (U M (r,t)- < U M (r,t) >f > 

2 



(i^)VI <fewi)in 



(21) 

where the correlation volume v o ~ A 3 and the integration over the variables y' and z' was 
restricted by the transverse size of the device, L±_ which is typically about few mm. The 
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contribution of correlations is about (A/r s ) 3 and remains much less than unity for all particle 
densities considered, therefore it does not affect the final estimation in Eq. (12 ip. 

For very low concentrations with r s ^> L, assuming constant p s (x) ~ r~ 3 all over the 
device one finds energy variation of the order of 

< (U M (r,t)- < U M (v,t) >) 2 >^ (j^—) ^ln^, (22) 

which cannot match the second term of Eq. ( fl3|) squared at any reasonable value of x m < L. 

This means domination of SP contribution for these concentrations in accordance with the 

mean-force analysis in the previous section. 

For intermediate concentrations such that It <r s -C L the equilibrium particle density 

may be used. Substituting Eq. (Tl6l) into Eq. (121]) and comparing the result with the second 

term of Eq. ( 1T31) squared one comes to the criterion 

f I \ 1/4 1 
x m ~ r s [ - ~ (0.05-0. l)r 8 (23) 

with the coefficient for r s by factor 0.25 less than that in the mean-force criterion ( fl8l) . 
When x m is less than the value given by the criterion ( 1231) the SP contribution dominates 
the injection process, in the opposite case the SP term becomes negligible comparing to the 
stochastic contribution to the energy so that the MF mechanism of injection prevails. 

For high particle concentrations with r s < If, assuming again constant p s (x) ~ Z^ 3 the 
energy variation of the form ( 122]) obtains. In the considered parameter range it dominates 
strongly over the SP term for all reasonable values of x m so that the MF contribution prevails 
in accordance with the mean-force analysis in the previous section. 

Summarizing analysis using both the mean-force and energy fluctuation criteria it is 
apparent that the SP approach to the injection is in any case relevant for low concentrations 
of injected particles p s < L~ 3 while the MF approach is valid for high concentrations p s > Z^ 3 
whereas individual contributions of single particles and their images are lost in the collective 
mean field. In the intermediate range L~ 3 < p s < Z^ 3 the two mentioned criteria bring about 
somewhat different estimations of the crossover from the SP to the MF domination regime, 
Eqs. ( fl8|) and ( |23|) . that can be roughly compromised by the criterion 

x m ^ 0.2r s , (24) 

which will be used in the following analysis. The above specified microscopic conditions im- 
ply restrictions on the macroscopic variables and system parameters which define description 
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Electric field [V/m] 



FIG. 1: Position of the potential maximum of the Schottky barrier, ^rn? 

and 0.2 of the mean 

distance between injected charge carriers in the MF approach, r s , as functions of the electric field 
at the interface for A = 0.225 eV and T = 300 K. 

of injection in terms of the SP or the MF approach when x m is smaller or larger than the 
value given by Eq. ( 124"1) . respectively. The corresponding macroscopic conditions for fields 
and voltages are revealed in the next section. 

D. Macroscopic criterion for the validity of the mean-field and the single-particle 
concepts 

Let us consider the intermediate range of concentrations of injected charge carriers, 
Ps(0) < Substituting expressions for r s = Ps(O)^ 1 / 3 , Eq. ([7]), and for x m , Eq. ffTUj) . into 
the criterion (1241 we arrive at a transcendent equation for the field F s (0) 



A graphical solution of the above equation for typical parameters of ITO and organic semi- 
conductors at room temperature (see Table [I]) is shown in Fig. [1] which exhibits two points 
of intersections. The dependence of the two corresponding field magnitudes on the bar- 
rier height A is shown in Fig. [2J The MF approximation is expected to be valid in the 
range of interface fields smaller than F m i n or larger than F max , where 0.2r s < x m , with the 
SP approximation prevailing between these lines as indicated in Fig. [2] with account of the 
assumed restrictions It < r s < L. 




(25) 
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FIG. 2: (Color online) Two field solutions of the equation x m = 0.2r s are shown with solid lines 
which delineate the areas where the MF and the SP concepts are valid. The dashed line Fu m 
restricts from above the field values calculated under the assumption of Boltzmann statistics. At 
the solid lines Ft and Fl the equations r s = It and r s = L are, respectively, satisfied. Thereby the 
area above Ft is described by the MF approximation. In the hatched area below Fl the continuous 
description of the transport and injection fails. To the left of the dash-dotted line at A = 0.16 eV 
the inequality 0.2r s < x m holds always. The other dash-dotted line at A = 0.38 eV marks the 
intersection of the lines F m i n and Fl. 

The area of high carrier densities p s (0) > l^ 3 is indicated by the line Ft where the 
condition r s = It is fulfilled. Above this line determined by temperature the density of 
injected particles is so high that individual image forces are strongly overwhelmed by the 
mean field. Note also the line Fi im indicating the border at which the Boltzmann statistics 
used in the MF consideration is violated. 

The preceding analysis in this section implied treatment in terms of the continuous par- 
ticle density and mean field. For the low carrier densities p s (0) < L~ 3 these quantities 
make no sense. The line Fl indicates the border where the condition r s = L is fulfilled. 
Below this line determined by the device thickness the injected particles are so sparse that 
the continuous description fails but the SP approach applies. This hatched area in Fig. [2] 
overlaps with the above mentioned SP area but has another physical meaning. Here the 
Richardson-Schottky formula^ for the SP thermionic emission should apply or more sophis- 
ticated concepts adapted to the hopping conductivity typical of organic systems^^ 1 ^— be 
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implemented. In the upper SP region above the lines F min and Fl the continuous description 
still works but the injection is dominated by the SP mechanism. 

To the left of the characteristic value of the injection barrier of A ~ 0.16 eV, where the 
lines F min and F max merge, the MF approach is valid for all fields between the values Fl 
and Fu m . To the right of the characteristic value of the barrier of A ~ 0.38 eV, at which 
the lines F min and Fl intersect, the SP approach is valid for all fields up to the value of 
Ft- In the region of barrier values 0.16 eV < A < 0.38 eV the chart becomes nontrivial. 
Namely, by increasing the electric field from the very low values, firstly the only SP approach 
is applicable, then MF approach becomes applicable too and describes injection, later the 
reentrant SP regime prevails due to domination of individual image forces, and by further 
field increasing the reentrant MF regime is realized at high carrier concentrations when 
reaching the lower one of the fields Ft or F max . The drift- diffusion equation (j3J) applies 
thereby for the description of charge transport inside the organic semiconductor in the 
whole region F L < F s < F Um . 

Note that all the characteristic field lines presented in the chart on Fig. [2] result from 
estimations like the criterion ( 12^1) and, hence, are not sharp boundaries but rather indicate 
the location of transition regions between the SP and MF domains. From the derivation 
of the criteria ({TBI) . (1231) and (pMI) it is clear that the uncertainty of the numerical factor 
in Eq. (125]) may cause only a week logarithmic change in the characteristic fields. The 
consequent uncertainty of the positions of the above field lines characterizes the width of 
the transition regions which appears to be much less on the logarithmic scale than the SP 
and MF regions themselves. 

Now, using the obtained threshold values of the fields in Fig. [2, we have to establish 
the range of voltages and barrier heights where the relevant field strengths can be achieved. 
Note that changing the injection barrier A for a fixed electrode means considering organic 
materials with different positions of the HOMO-level. 

Combining the solution (J5j) with the boundary condition (JTj) one finds that the field at 
the interface is negative in equilibrium (j = 0) or when the applied voltage is small*^ 1 ^. 
The reason for this is that the field has to compensate the positive diffusion current caused 
by the huge difference in charge-carrier concentrations between the different sides of the 
interface. When a positive bias is applied, the field changes its sign at some point inside the 
dielectric usually called a "virtual electrode"-. In the parameter range where the electric 
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field is negative at the interface the Schottky barrier concept does not apply because a single 
charge carrier cannot escape from the potential well. In this case, however, the collective 
injection can occur which thus could be described within the MF approach. The boundary 
of this parameter region is delineated by the condition that the field at the interface equals 
zero, or, in other words, that the virtual electrode coincides with the physical interface. For 
a given value of the barrier the current magnitude at which this occurs may be found exactly 
using the solution (jSJ) and reads 
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where zq ~ —1.02 is the first zero of the Airy function Ai'(z) (Ref.— ). The current- voltage 
(I-V) relation for a single ITO/organic interface is^ 
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where the constant F s q(0) has to be determined from the boundary condition (0) at the in- 
terface in the equilibrium case j = while the constants F s (0) and C\ have to be determined 
satisfying Eq. ((7j) with the solution ([5]). Using the current value jo and the current- voltage 
relation (127")) a magnitude of the voltage Vo(A) may be obtained which defines the line 
F s (0) = on the V — A plane (see Fig. [3]). Using the characteristic field values of the 
Fig. |2j also the voltages V min , V max , Vu m , Vt and Vi which separate the MF and SP regions 
can be found by adjusting the current until the respective field values F min , F max , Fu m , Ft 
and Fl are achieved at the interface as is presented in Fig. [31 The above calculations using 
continuous particle densities are applicable as long as the inequality r s < L roughly holds 
that is above the line Vl below which only SP consideration makes sense. The line Vt 
crosses the line V max so that the SP dominated region is bound from above by the volt- 
age min {V max , Vt}- Above this voltage the MF approach is valid until the voltage Vu m is 
achieved at which it fails together with the Boltzmann statistics. 

Though it looks counterintuitive, the continuous MF description of injection may be 
valid for low voltages while the discrete SP description, expected to work at low injection 
rates, applies also for rather high voltages. The reason can be seen in the very small escape 
distances x m at high voltages for which x m r s , if the injection barrier A is not low, so 
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FIG. 3: (Color online) Two solid lines, V m in and V max , outline areas of validity of the SP and the 
MF approaches. The increasing dashed line Vt further bounds the SP region from above. The line 
Vn m bounds from above the region where the Boltzmann statistics applies. The decreasing dashed 
line Vo, shows the voltage at which the field at the interface vanishes. The line Vl delineates the 
boundary of the hatched area where the mean field treatment loses any sense. 

that the individual image force on the charge carrier dynamics dominates. At low voltages 
and low barriers, however, the density of injected particles is diffusion- driven and therefore 
relatively high which reduces the interparticle distance r s to values small compared to the 
distance for the individual particle escape x m and entails domination of the mean field. 

III. MODIFIED MEAN-FIELD BOUNDARY CONDITIONS FOR A SINGLE IN- 
TERFACE 

The deficiency of the MF approach to the description of the injection within the SP 
parameter range in Fig. EHbelow the lines min {V max , Vt} and above the lines max {V min , Vl} 
is the missed contribution from the individual image forces of each single charge carrier. This 
strong but short-range deviation from the mean field near the interface may be accounted 
for as a dipole layer of a characteristic width x m as long as the latter distance is much 
less than r s , which is true in the mentioned region. Indeed, this layer contributes to the 
shift in energy of each injected charge carrier through the work performed by the individual 
image force, which results in the shift in the electrostatic potential at the interface given by 



18 



8(f) schi Eq. (ITT]) . Accounting for that in the mentioned SP region the boundary condition 
(requirement 1. of continuity of the electrostatic potential in the Section III A I) has to be 
replaced by the condition 

- C (-O) = 8(j) Sch . (28) 

Applying this boundary condition at x = the dipole layer is formally considered as an 
infinitesimal sheet so that the drift-diffusion equation fl3j) applies for x > 0. 
Taking this into account the modified boundary condition reads 



Ps(0) 

— = exp 



kT + kT e r [J) 



1_ e*F s (0) 
0.2r s y kT\] A7T6 e s 

in the SP region restricted by the inequality max{V m i n , Vl} < V < mm{V max , Vr} with 
b = max{i m , 0.2Zj>}. An interpolation factor (1 — 6/0. 2r s ) introduced in the last term 
provides smooth switching of the dipole layer contribution at the boundary of the SP region 
which ensures that the position of this boundary is defined consistently when evaluated from 
inside and outside of this region. Indeed, crossing the line V m i n or V max from the SP region 
to the MF region the length x m exceeds 0.2r s and can no more be considered as the dipole 
layer thickness. In this case, the image force averaged over all injected carriers and their 
images is accounted for in the mean field value F s (0). The same occurs also at the line 
Ft according to analysis in Section III CI Thus, the last term in the brackets in Eq. ( 129]) 
following from the potential difference (128]) . has to vanish at the boundary of SP region 
which is provided by the above introduced interpolation factor. 

Now Fig. [3] has to be reconsidered on the basis of the modified boundary condition (|29j) . 
It is obvious that the lines Vo, V m i n , V max and Vr do not change their positions since, at these 
lines, the boundary condition (J7]) still holds. The equation r s = L, however, is affected by 
the modified density of injected particles in Eq. ( 129]) . This is followed by the revised value 
of the field Fl producing the new borderline Vl- The MF approach applies now also in the 
SP region above the line Vl, however, with account of the dipole layer in Eq. ( 129]) . therefore 
it will be called from now on a modified mean-field (MMF) region. The genuine SP domain 
where the SP concept only applies shrinks now to the region below the Vl line. 

Finally, the generalized boundary conditions (J7|) and (129]) apply on the whole V — A plane 
below the V[ im and above the Vl lines and can be used together with standard continuous 
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FIG. 4: (Color online) Solid voltage lines V m i n ,V max ,VT and, partly, Vl outline the MMF area 
where the MF concept with the boundary condition f|29|) applies in contrast to the MF area where 
the boundary condition ([7]) holds above and below the line Vq. The latter line shows the voltage 
at which the field at the interface vanishes. The MF approach fails completely only in the hatched 
SP region under the line Vl- The line Vu m bounds from above the region where the Boltzmann 
statistics applies. 

equations for the field and charge densities while the RS or FN single-particle concepts are 
expected to hold below the line Vl- To make these results comparable with experiment we 
address below two examples of practically relevant two-electrode devices. 

IV. MODIFIED MEAN-FIELD APPROACH FOR TWO-ELECTRODE DEVICES 

In the case of the organic layer sandwiched between two electrodes, located at the posi- 
tions x = ±L/2, multiple images of an injected hole will contribute to the potential in the 
SP approach. Consequently, the second term in Eq. (jH]) has to be replaced by the energy of 
charge interaction with both electrodes derived in the Appendix. We have to account now 
for possible injection of charge carriers from both electrodes. If the mean field F s becomes 
positive at the left interface at x — —L/2 or negative at the right interface x = L/2 (which 
never occurs at the same time^i) the SP Schottky mechanism of injection may become rele- 
vant at the respective electrodes. In the spirit of our modification of the boundary condition 
at the single interface x = in Eq. (I2"gj) - (|2"§|) . the boundary conditions at x — ±L/2 should 
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be modified by the potential discontinuities 5(f)^ ch if the distance from the potential max- 
imum to the respective electrode L/2 — \x m \ < 0.2rf with rf = [p s (±L/2)]~ 1 ^ 3 being the 
mean distance between injected particles near the electrodes at ±L/2. In this case, the 
mean field may be assumed constant within the distances rf from the respective electrodes 
and the following approximations used for the energy profile near the electrodes 

U± (x) = A ± - eF s (±L/2) (x T L/2) 



e 2 



+ iroK2 + zj + n2-ij +27 J' (30) 

where A ± denote the injection barriers at x — ±L/2, 7 is Euler's constant and ip(x) is 
the digamma function^. The expression in square brackets reproduces the behavior of the 
Coulomb potential in Eq. fl9]) when approaching one of the electrodes at x = ±L/2 but is 
symmetric and has a maximum at x — 0. Note that the criterion (|24p derived using the 
energy (Q for a single interface accounts, in fact, for the second electrode and can therefore 
be applied to the two-electrode devices. 

According to the approximations ( 1301) the position of the maximum is defined by one of 
the equations 

* G + T) " * S " T) = ^f^F. (±i/2) , (31) 
taking the upper sign, if F s (L/2) < 0, or the lower sign, if F s (—L/2) > 0. The respective 
potential discontinuities are given then by the formula 

e6^ ch = A ± -U ± (x m ) , (32) 

which determines them as functions of the field at the respective interfaces. 

The boundary conditions 1.-3. of Section Hi Al modified with the potential discontinuities 
(12 8p and ( 132]) at both interfaces can now be written as 



Ps(±L/2) 
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exp 



A ± eF s (±L/2)l± F e s 



. i b ± \ e5(j)± ch 



kT 1 kT e± 

0.2rfJ kT J' (33) 
where quantities identified with superscript ± denote here and below the parameters of the 
two electrodes contacted at x — ±L/2. Note that the last term in the exponent appears only 
inside the MMF area on the V-A plane whose boundaries are determined self-consistently 
by the criterion L/2 — \x m \ = 0.2rf using the formula (1331) itself. The cases of symmetric 
and asymmetric devices are considered separately below. 
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A. Modified validity chart and current-voltage characteristics of a unipolar sym- 
metric organic device 

We examine first the case of a symmetric unipolar device consisting of two ITO-electrodes 
and an organic layer in between providing p-type conductivity. At zero bias the diffusion- 
mediated electric field is directed outwards at both electrodes, i.e. is negative at x — — L/2 
and positive at x = L/2. In this case the Schottky-type contribution in the formula (133]) 
vanishes and the MF boundary conditions of the type of Eq. ([7]) apply at both interfaces 
as long as the barrier is not too large which is in agreement with the chart for the single 
ITO/organic interface (Fig. H]). Note, further, that for a positive bias, a positive field always 
prevails at the right, collecting electrode so that the Schottky-type barrier lowering does 
not arise here at all and this electrode remains in the MF regime as long as the continuous 
approach is valid. 

On the other hand, when the bias is high enough, the field at the left, injecting electrode 
becomes positive and may exceed the characteristic field F m i n in which case the modified 
boundary condition ( )33l) includes the Schottky lowering term. Since S(pg ch is not known 
as explicit function of the field F s (—L/2) the equation (133]) should be solved together with 
Eqs. f l30|) -f l32|) . Then the mean interparticle distance r~ can be evaluated and the equation 
L/2 — \x m \ = 0.2r~ solved with respect to the field F s (—L/2) resulting in the positions of 
the lines F min and F max on the F — A plane for the injecting electrode. Additionally, the 
characteristic field Fu m can also be found equating the exponent in Eq. ( )33l) to zero, and the 
field Ft can be determined from the equation r~ = It- Determination of the line Fl, where 
the continuous approach fails, is a bit less straightforward. When the charge carriers are so 
sparse in the device that the mean distance between them equals the distance between the 
electrodes their concentration is defined by injection from both electrodes. Therefore the 

- - ——1/3 

mean particle density p s and the corresponding mean distance r s = p s should be used for 
the definition of Fl from equation f s = L. Integrating the Gauss equation over the device 
thickness one finds that 

p s = ^ [F s (L/2) - F s {-L/2)\ . (34) 

Thus the procedure of Fl determination is as follows. Below the line F min the MF boundary 
conditions, Eq. ([7]), are satisfied at both electrodes with Eq. ([5]) and the current is increased 
until equation f s = L is fulfilled. Above the line F m i n the boundary condition ( 1331) should 
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FIG. 5: (Color online) Solid lines V m i n , V max , Vt, and, partly, Vi outline the MMF area for a sym- 
metric ITO/organic/ITO device where the SP contribution to the MF approach must be accounted 
for. The dashed line indicates the voltages Vo at which the field at the left interface vanishes. The 
line Vii m bounds from above the area where the Boltzmann statistics applies and the line Vl the 
hatched area where the continuous description fails completely. 




be satisfied at the left electrode and the MF boundary condition (J7|) at the right electrode 
while the current is adjusted until the requirement f s = L is fulfilled. The mentioned lines 
delineate together the MMF region on the field-barrier plane similar to Fig. [2j This chart, 
however, plays only an auxiliary role and will not be presented here. 

Having determined F min , F max , Fn m , Ft and Ft as functions of A~ one can calculate, 
using the solution (j^D, the respective voltages V , V min , V max , Vu m , Vt and Vp which delineate 
the MMF region. To this end, the current is varied as parameter until the respective field 
magnitude is achieved at x — —L/2; the voltage then follows by direct integration of the 
field over the device thickness^i^. Thereby the free constants C\ and Ci in Eq. ([5]) are 
obtained by satisfying the conditions (|33|) and ([7]) at x = —L/2 and x = L/2, respectively. 
This method, however, appears to be numerically unstable because of fast oscillations of the 
Airy functions providing multiple solutions for the constants. Alternatively, Eq. (j3J) can be 
directly numerically integrated using the same boundary conditions which proves to be a 
robust procedure. 

The resulting chart on the V — A~ plane is depicted in Fig. [51 Apparently, this chart 
is very similar to that of the single interface depicted in Fig. 0] except the voltage region 
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FIG. 6: I-V characteristics for symmetric organic devices with and without account of charge 
carrier discreteness are shown by solid and dashed lines, respectively. The barrier magnitude is 
changed from the upper to the lower curve as is shown in the plot. For the barrier values below 
0.2 eV and smaller the solid and the dashed lines cannot be distinguished. 

below 0.1 V where the Vq and V m i n curves go down a bit steeper. This means, on the one 
hand, that the injecting interface dominates the behavior of a symmetric device. On the 
other hand, this similarity shows that accounting for the collecting electrode in the energy 
profile of a single particle, Eq. fl30|) . does not add much to the simpler formula for a single 
interface, Eq. ([9}. 

I-V characteristics for different injecting barrier values are plotted in Fig.[6]for a symmet- 
ric device, employing typical parameters of ITO and organic semiconductors (see Tabled]). 
The essential contribution of the SP Schottky barrier lowering in the wide region of voltages 
is evident. 

B. Modified validity chart and current-voltage characteristics of an unipolar asym- 
metric organic device 

Asymmetric devices are considerably different from symmetric ones due to the presence 
of the built-in potential, Vu, which is defined by the difference between the work func- 
tions of the two electrodes, eVu = E\ — E^ = E£ + k+ — E^ — k^. At equilibrium, 
due to equalization of the electrochemical potential in the entire system, an internal elec- 
tric field emerges. For that reason the electric field at one of the electrodes may change 
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its sign providing favorable conditions for the SP injection scenario already at zero bias. 
To account for this effect the field-barrier charts should be considered separately for both 
electrodes by solving equations L/2 — \x m \ = 0.2rf which result in corresponding charac- 
teristic fields for both electrodes. Consequently, Eq. should be solved implementing the 
boundary conditions self-consistently at x = ±L/2. Thereby potential discontinuities 
Wsch from E( l- © should be applied in Eq. §3§ as long as max{|F^ n (A ± )|, \Ff (A ± )|} < 
|F a (±L/2)| < min{|F± a:r (A ± )|, |F± m (A ± )|, |F±(A ± )|}. Note that the fields F± are defined 
simultaneously from the mean value of the carrier density using Eq. as for the symmet- 
ric device. As soon as the characteristic fields are known the corresponding characteristic 
voltages may be evaluated as was described in the preceding section. The line Vl is obviously 
common for the V — A 1 * 1 charts of the left and right electrodes. 

Application of the self-consistent boundary conditions f )33|) may look a bit cumbersome 
but it is substantially simplified by the fact that the sole maximum of the energy (130|) . x m , 
occurs at each moment only near one of the two electrodes. Thus, if the external bias has 
the same sign as the built-in potential then the SP contribution may arise at the injecting 
electrode only. If the external bias has the sign opposite to that of the built-in potential the 
SP contribution may arise first in equilibrium at the collecting electrode due to the internal 
electric field while the injecting electrode is described in the MF approach. By the increasing 
external field the collecting electrode goes over to the MF regime too and by further field 
increasing the injecting electrode goes over to the SP regime where the discontinuity 54>s c h 
applies. 

To investigate possible effects induced by the presence of the built-in potential, the case 
of an asymmetric unipolar device consisted of the organic layer sandwiched between ITO 
and Al electrodes is considered. For the ITO electrode at x — —L/2 and Al electrode at 
x = L/2 the built-in potential equals -0.4 V (see the electrode parameters in Table HT|) . The 
V — A - chart for this system is presented in Fig. [7] where the behavior of the right electrode 
is also comprised. Changing the injection barrier A - for fixed electrodes means considering 
organic materials with different positions of the HOMO-level. This implies, in turn, the 
change of the injection barrier at the other electrode, A + , to the same extent. This allows 
one to show the area where the SP injection mechanism is relevant at the right electrode on 
the same V — A~ chart. 

The Vq curve for the asymmetric case, where the field vanishes at the left, injecting 
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TABLE II: Parameter values for collecting electrodes (see Ref. 
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FIG. 7: (Color online) Solid lines V m i n ,V max ,VT and, partly, Vl outline the MMF area for an 
asymmetric ITO /organic/ Al device where the SP contribution to the MF approach must be ac- 
counted for at the left, injecting electrode. The upper dashed line indicates the voltages Vq at 
which the field at the left interface vanishes. The lower dashed line delimits the low- voltage region, 
where the SP contribution to the MF approach must be accounted for at the right, collecting Al 
electrode. In the hatched SP region bound by the line Vl the description of the device with the 
continuous carrier density fails completely. The Boltzmann statistics is valid below the line Vu m . 




interface, nearly coincides with that of the symmetric case (cf. Fig. \5§ for barrier values 
below A~ = 0.22 eV. Above the value of 0.35 eV it coincides virtually with the lower 
boundary of the MMF region given by the V m in curve because F min value becomes equal 
to zero. When A~ > 0.25 eV, the Vq curve differs substantially from the corresponding 
line for the symmetric device and saturates at about the value of —Vu- The Vu m curve 
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FIG. 8: (Color online) I-V characteristics for asymmetric organic devices with and without account 
of charge carrier discreteness are shown by solid and dashed lines, respectively. The barrier mag- 
nitude is changed from the upper to the lower curve as is indicated in the plot. The solid circle 
on the characteristic for A = 0.4 eV marks where it enters to the low- voltage region (dash-dotted 
line) where the mean-field approach is no more valid. 

which delineates the upper boundary of the chart almost coincides for the symmetric and 
asymmetric cases. Thus, the presence of the built-in potential reveals itself in the extension 
of the low-voltage area, where the Schottky correction can be neglected at the injecting 
electrode, up to approximately — V hi . On the other hand, this correction should be included 
at the collecting electrode at almost all voltages below — V u as is indicated by the lower 
dashed line. 

I-V characteristics for different injecting barrier values are plotted in Fig. [HJfor the asym- 
metric device, employing typical parameters of ITO and organic semiconductors 32 . As in 
the symmetric case, the Schottky lowering of the barrier enhances the current by orders of 
the magnitude for high voltages compared to the pure MF description^ 2 .. This results in 
the disappearing of the second exponential increase in the case of asymmetric devices and, 
generally, in an earlier transition to the space charge limited regime in both symmetric and 
asymmetric cases. Note that account for the charge carrier discreteness does not change the 
I-V curves below — Vu in Fig. [HJ 
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FIG. 9: (Color online) Measured I-V characteristic of ITO/OCiCio-PPV/Au structure (squared 
curve) and its approximations within different theoretical models. The best fitting parameters in 
the MMF approach are indicated in the plot. 



V. COMPARISON WITH EXPERIMENT 

To evaluate the presented model, the results for an asymmetric device are compared with 
experimental data obtained on a unipolar device consisting of a single poly[j>-(2-methoxy- 
5-(3,7-dimethyloctyloxy)phenylene-vinylene] (OCiCio-PPV) layer of thickness L = 100 nm 
sandwiched between ITO and Au electrodes. The characteristic energy of the HOMO level 
of the employed OCiCio-PPV is 5 eV and the hole mobility fi s is reported between 5 x 10~ 7 
and 3 x 10~ 5 cm 2 /(V s) or even larger— ~—. Note that for an organic semiconductor the 
mobility is typically scattered by orders of the magnitude^. It is known also that the work 
function of ITO is sensitive to the cleaning procedure and thus, it may vary from 4.7 to 5 eV 
(up to the energy of the HOMO level). Additionally, a dipole layer may possibly emerge 
at the metal/organic interface, typical for these systems^ 1 ^. Therefore, we consider the 
injection barriers A~ and A + as independent fitting parameters. 

In Fig. |9]the measured I-V characteristic of the ITO/OCiCio-PPV/Au structure is shown 
as well as the I-V dependencies calculated for different model approximations. One can see, 
that the advanced MMF approach gives a good approximation for the measured I-V curve 
in the wide range of applied voltages with the best fitting parameters indicated in the figure. 
At the same time, the dashed curve for the MF approach with the same fitting parameters 
underestimates the current by orders of the magnitude at high voltages. Note that the 
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widely used Crowell and Sze model 2 -^ (the dash-dotted curve) fails to describe the I-V 
characteristic at voltages below —Vu- 

Notice, however, that the best fit was found for slightly higher magnitudes of the injection 
barriers than the largest values reported in the literature: A~ ax = 0.3 eV and A+ ax = 0.7 eV 
(the work function of Au is 4.3 eV, see Table HI])- This deviation of the barrier values may 
originate from the uncertain position of the HOMO level in OCiCio-PPV depending on 
details of the material production as well as from the above mentioned possible dipole layers 
at the electrode/OCiCio-PPV interfaces. Also, the exact shape of the I-V curve cannot be 
fully approximated because of the fact that the proposed model still misses some important 
features of insulators, in particular of organic semiconductors. Those are the realistic DOS 
shape including energetically distributed trap states^ 1 ^ and the field- and concentration- 
dependent carrier mobility—^ characteristic of organic semiconductors. Remarkable devi- 
ation of the theory from the experiment below the current densities of 10~ 4 mA/cm 2 may 
have origin in some other physical mechanisms as for example leakage effect or impurity 
doping. Nonetheless, a sufficiently good matching of experiment and MMF model is found 
over 2 orders of magnitude in the applied voltage and 7 orders of magnitude in the current 
density. 



VI. DISCUSSION AND CONCLUSIONS 



In this work, we have advanced boundary conditions for a continuous description of charge 
injection and transport in dielectrics in terms of mean electric field and charge density ac- 
counting for discreteness of charge carriers. The latter becomes important for a wide range 
of values of injection barriers and voltages applied where individual image force dominates 
the injection process. Regions of parameters, where charge carrier discreteness plays a deci- 
sive role, are approximately delineated in the V — A charts for symmetric and asymmetric 
conductor /insulator /conductor devices using an organic semiconductor as example for an 
insulator sandwiched between two electrodes. Implementing generalized boundary condi- 
tions allows for the description of charge transport in such systems in both low injection 
and high injection regimes including space charge limited transport. The advanced theory 
applies to a wide class of devices where injection of charge carriers into dielectrics takes 
place, particularly, to organic light-emitting diodes, thin film ferroelectric capacitors, etc. 
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Though in this work only the problems of simple plain geometry were treated, the developed 
approach is not intrinsically one-dimensional and can be applied as well to complicated two- 
and three-dimensional configurations of conductor /insulator interfaces as, for example, those 
of organic field-effect transistors, however, with duly adapted boundary conditions fl3"3"|) . 

The presented description of injection applies directly to inorganic crystal insulators and 
wide-gap non-degenerate semiconductors as well as to the very narrow-band insulators and 
semiconductors as was indicated in Section III Al When applying this theory to the wide- 
gap organic semiconductors one should take into account specific features of these media. 
Being disordered semiconductors with typically hopping conductivity, organic systems differ 
from inorganic crystals mostly in two respects concerning the injection process, namely, 
they 1) possess as a rule a relatively wide (Gauss) DOS distribution without sharp band 
edges^^^i^ 1 ^ 1 ^ 1 ^ and 2) possibly exhibit long-range correlations of the random energy 
levels^. 

Considering point 1), the absence of a sharp band edge makes the injection barrier an ill- 
defined parameter. If the barrier A is defined as a difference between the Fermi level in the 
injecting electrode and the center of the density of states distribution at the HOMO-level as 
it was done in Refs.— iI2i2L2& the injection into the "gap" states below the barrier is thereby 
implied because of numerous states available at least to the depth of A — cr, where a is the 
width of the Gauss DOS. Nevertheless, a nontrivial fact is the possible injection into the tail 
states well below the barrier— (to the depth of about A — a 2 /kT) which results in 
the field-dependent mobility of carriers^ or in the allegedly possible space charge regime^ 
for high injection barriers. The problem of the concept of injection into the deep tail states 
is that the injection is considered as a single-particle process ignoring the fact that these 
states may be occupied. Particularly problematic becomes, therefore, a combination of the 
SP injection with the drift- diffusion equation^ 8 at high injection levels because in the space- 
charge regime the density of injected charge is especially high near the injecting electrode. 
On the other hand, the concept of injection into the tail states may be well applicable at 
low injection levels, particularly, in the genuine SP region of the V-A chart of Fig. [7J In our 
simulation of the experiment in Fig.[9]we have assumed a narrow-band approximation for the 
HOMO-level ignoring the finite width of the DOS and, thus, the disorder effect on injection. 
This might explain the deviation of the theoretical characteristic from the experimental 
one at low voltages. In fact, to ignore the role of disorder on injection the inequality 
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a <C kT should be satisfied which is not typically the case for organic semiconductors. 
Proper integration of the finite DOS width in our self-consistent concept is currently in 
progress. 

A rather elaborated many-particle description of charge injection and transport in organic 
semiconductors with uncorrelated Gauss disorder was developed recently in the Ref.— where 
boundary conditions similar to our MMF ones were used to account for the image charge 
effect. Two models were presented and compared in this paper, a discrete model of 3D 
hopping over the sites of a cubic lattice and a continuous ID description in the spirit of the 
drift-diffusion concept, which exhibit good agreement with each other. However, in both 
models the individual image effect was substantially overestimated. In the 3D approach the 
individual image contribution was accounted for at all internal sites of the lattice in addition 
to the mean field derived from the charge averaged over the planes parallel to electrodes. As 
was discussed above in Sections [III and [Till the individual image effect reveals itself as strong 
deviation from the mean field only in the vicinity of the electrodes due to its short-range 
nature. In fact, close to the electrode surface the individual image force may dominate over 
the mean field. In contrast to this, far from the electrode a charge carrier "observes" not 
only its own image but all the other injected particles with their images which together 
constitute the mean field. In the latter case there is no reason to single out the carrier's own 
image which would mean duplication of the image account. This problem, was discussed in 
detail in Refs.— ~— where such duplication was explicitely excluded. In the continuous ID 
model of Ref. 37 , the line Vo was chosen as a boundary between the MF and MMF region, 
where the barrier lowering due to the image effect is included. This line is close to our 
lower boundary of the MMF region V m i n on the V- A charts in Figs. HJ [5] and [7] for a certain 
range of barrier values. This approach, however, also overestimates the image effect for low 
barriers and misses the fact that below some threshold barrier value of about 0.2 eV the MF 
approach is always valid independently of the field direction at the electrode. 

The property 2) of organic semiconductors which may be relevant to statistical grounds of 
our concept of injection, namely, the correlated disorder of the energy levels at different sites 
is not so well established and is still disputed in the literature. The idea of the correlated 
disorder was suggested in Refs.— to explain the field dependence of mobility In /x ~ y~F 
in a wide region of field values. We note, firstly, that this single-particle concept is based 
on the arbitrary hypothesis of independently and randomly oriented dipoles at each lattice 
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site^. Secondly, such field dependence of the mobility was also explained within another 
single-particle model considering the injection in the deep tail states of the uncorrelated 
Gauss disordered energy levels^-. Many-particle models assuming uncorrelated disorder 
combined with concentration dependence of the mobility^^i can provide good agreement 
with experiments, too. Pasveer et al. concluded that there is no necessity to account for 
correlated disorder at least around room temperature^. Hence, possible correlations in 
charge positions due to the correlated disorder are not regarded in our study. 

Considering the Schottky-barrier lowering we have to note the general problem of this 
concept caused by the too large distances between the electrode and the potential barrier 
maximum at low voltages. In fact, this value is restricted for different physical reasons 
depending on the material type. In crystalline dielectrics this concept, implying usually 
ballistic overcoming the potential barrier by a charge carrier—, makes sense as long as 
the mentioned distance does not exceed the mean free path which may be large enough. 
On the other hand, in organic semiconductors characterized by hopping conductivity the 
ballistic description fails because of typical hopping distance of 0.1—1 nm. In this case, 
more sophisticated descriptions of injection like multiple random hopping of a single charge 
carrier in the SP potential may be relevanlr^ 1 ^. Applicability of the latter model seems 
to be restricted to the SP region in the charts of Figs. HJ and [7] where the effect of the 
other carriers is reduced. After all, we would like to stress that in the formulation of the 
thermodynamic boundary conditions 1.-3. of Section Hi Al together with Eq. (J28l) the ballistic 
mechanism of charge injection was not assumed. 

Finally, we considered in this paper the corrections to the self-consistent mean-field de- 
scription of injection due to the effect of discreteness of charge carriers. It is apparent, 
however, that below the line Vl, where continuous description is no more valid, the pure SP 
approach will also be distorted because of the long-range interaction with the other injected 
charge carriers at the distance L or larger. This is a sort of MF corrections to the SP picture 
which is yet to be elaborated. 
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Appendix: Energy of a point charge between two plane electrodes 



When a point charge q is located in the dielectric medium with relative permittivity e s 
between two plane electrodes taking positions at x = ±L/2, the Coulomb potential used in 
the formula ([9]) for the single interface does not satisfy anymore the boundary conditions 
for the electrostatic potential at the electrodes: ip(x = ±L/2) = (no voltage applied). To 
satisfy these boundary conditions the method of images may be used~ which results, in the 
case of two electrodes^, in the potential expression 

q ^ (-1)" 



ip{x,y,z) 



(A.l) 



4vre e s ^ - Xn ) 2 + y 2 + z 2 ' 

where (xo,0, 0) is the position of the point charge and x n = xq(— l) n + nL with n = 
±1, ±2, . . . are the x— coordinates of its images. 

The force, exerted upon the point charge, results from the interaction with image charges 
to the left and to the right of its position: 



f(xo) 



Q 
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2q 2 x 
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oo 
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where a summation formula from Ref. 
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1 Xo 

2 L 



Xq 

L 



(A.2) 



was used and ip(x) denotes the digamma function 42 . 
The energy of the point charge interaction with electrodes, U q , can be obtained by integra- 
tion of the force ( 1A.2I) considering the relation f(xo) = —U'Jxq). The integration constant 
U q (0) should be chosen so that, in the limit \x + L/2\ L, the energy U q (x ) reduces to 
the case of a single interface, Eq. ([9]). Thus, the expression 



U q (x ) 



167re e s L 

results, where 7 = 0.5772 is Euler's constant 42 . 



(A.3) 
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